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ABSTRACT 

Pneumatic  hammer  instability  of  an  inherently  compensated  thrust  gas  bear¬ 
ing  was  analyzed  theoretically.  Vohr's  experimental  correlation  between  the 
pressure  loss  coefficient  and  Reynold's  number  was  used  calculate  flow  through 
tne  restrictor  instead  of  using  cne  nozzle  equaffcr.r  The  time  d-^wndent 
Reynold's  equation  was  solved  by  a  perturbation  analysis  for  small  axial 
oscillation.  Based  on  the  perturbation  analysis,  dynamic  stiffness  ard  damp¬ 
ing  coefficient  were  calculated.  Utilizing  these  and  Pan's  stability  criteria 
(Ref.  6)  stability  maps  were  constructed. 


1.  INTRODUCTION 


Externally  pressurized  gas  bearings  have  been  used  in  many  engineering  devices. 

It  is  well  known  that  in  order  for  the  bearing  tw  have  relatively  large  load  capa¬ 
city  and  stiffness  it  is  desirable  to  have  recessed  pockets  immediately  after  the 
feeding  holes.  This  causes  the  externally  pressurized  gas  bearings  to  be  suscep¬ 
tible  to  pneumatic  hammer  instability.  Analytical  investigations  on  this  subject 
were  made  in  References  1,  2  and  3. 

In  conventional  analyses  of  externally  pressurized  bearings,  nozzle  equations 
are  used  in  calculating  the  flow  across  a  restrictor.  The  dynamic  pressure  head 
resulting  from  expansion  through  the  restrictor  is  assumed  to  be  completely  lost 
when  entering  the  fct aring  film.  This,  however,  is  not  true  as  reported  in  Refer¬ 
ences  4  and  3;  a  measurement  of  pressure  at  the  restrictor  exit  indicates  that 
there  is  considerable  pressure  recovery.  It  was  shown  that  the  pressure  loss 
coefficient  can  be  correlated  with  the  Reynolds'  number  (Ref.  4);  a  linear  rela¬ 
tionship  is  chosen  for  simplicity. 

A  simple  thrust  plate  with  a  feeding  hope  at  the  center  and  a  recessed  pocket 
is  to  be  analyzed  based  on  the  above  pressure  loss  coefficient  correlation  for  the 
restrictor  flow  and  the  Reynolda'  equation  for  the  bearing  film.  Perturbation  analysis 
for  small  oscillation  about  the  equilibrium  position  will  be  performed.  Eased  on  the 
perturbation  analysis,  dynamic  bearing  stiffness  and  damping  coefficient  are 
calculated.  Using  the  stability  analysis  of  Ref.  6,  stability  maps  are  constructed. 
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2.  ANALYSIS  -  SMALL  PERTURBATION 


The  configuration  of  an  inherently  compensated,  hydrostatic,  circular,  thrust 
bearing  is  schematically  shown  in  Figure  1,  Gas  at  supply  pressure  pg  is  led 
through  the  feeding  hole  with  diameter  dj,into  the  recessed  pocket  before  enter¬ 
ing  the  bearing  film.  For  a  circular  bearing  it  is  convenient  to  use  the  polar 
coordinates.  If  we  further  assume  circular  symmetry,  l.e.  no  misalignment,  then 
the  radial  coordinate,  r,  is  tne  only  space  variable  required  to  describe  the  flow 
and  the  pressure  uisti'iUuLiou.  In  oru^r  to  facilitate  »  dvnsmic  analysis  let  us 
allow  the  bearing  to  have  small  axial  vibrations  about  its  equilibrium  position  and 
express  the  bearing  film  thickness  as 

h  -  C  +  e  cos  t  (2.1) 

or  in  dimensionless  form 

h  •-  1  +  e  cos  t  (2.2) 


where 

h  -  h/C 
e  •  e/C 

C  *  equilibrium  film  thickness 
t  -  tit  ■  dimensionless  time 
to  ■  frequency  of  vibration 


(2.3) 


We  have  'ssumed  that  the  vibrations  are  purely  sinusoidal.  Note  that  t,  the 
normalized  amplitude  of  vibration,  is  a  small  number. 


The  well-known  time-dependent,  isothermal  Reynolds’  equation  is,  in  dimensionless 

form, 


I  J_ 

r  3r 


(2.4) 


where 


t/  K 

P/Pa 

12pu 


a 


R 

C 


squeeze  number 


(2.5) 
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The  boundary  conditions  are 
at  r  *  1  p  ■  1 

■  (2.6) 

dR  . 

at  r  "  2R  ’  P  "  PR/pa  "  PR 

While  Eq.  (2.4)  governs  the  pressure  distribution  in  the  bearing  film,  addi¬ 
tional  pressure-flow  relationships  across  the  inlet  restrictions  at  r  *  rf  and  at 
r  ■  rD  are  required  for  the  solution  of  the  problem.  In  the  literature  (Ref.  1*2,3) 
the  well-known  nozzle  formula  is  used  to  calculate  the  expansion  of  air  from  p 

a 

to  p_  and  from  p_  to  p_  (assume  that  the  pressure  is  uniform  in  the  recessed  pocket; 
see  Appendix  A) .  If  the  pressures  calculated  according  to  the  nozzle  formula  are 
accepted,  one  automatically  assumes  that  the  velocity  iiead  resulting  from  expansion 
through  the  nozzles  Is  completely  lost.  This  is  not  so  because  part  of  the  velocity 
head  is  recovered  as  indicated  by  references  4  and  5.  In  fact,  if  we  express  the 
pressure  drop  at  the  entrance  in  terms  of  the  velocity  head, 

K'  ‘t  <2-7) 

where 


K1  *•  loss  coefficient 
p  **  downstream  gas  density 
V  «*  downstream  average  velocity 


In  Ref.  4,  K’  is  correlated  experimentally  with  the  Reynolds'  number  (Re)  which  is 
reproduced  in  Figure  2.  But  for  all  practical  purposes,  a  linear  relationship  be 
tween  K'  and  Re  is  satisfactory, 

K'  *  K  Re  -  K  ^  (2.9) 


where 


K  ■  constant  «  0.66  x  10-d 
L  -  typical  length  ■»  film  thickness 


(2.10) 


ent 


(2.11) 
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Applying  the  above  equation  to  the  inlet  of  the  recessed  pocket  we  have 


Ps  ' 


V 


(h 


+  V 


(2.12) 


Inside  the  recessed  pocket  we  can  write  another  Reynolds'  equation  applicable 

there  and  solve  for  the  pressure  distribution.  However,  past  experience  indicates 

that  the  gradient  of  the  square  of  the  pressure  varies  inversely  as  the  cubic  of 

the  local  film  thickness  (see  Eq.  (A-5) .  In  most  applications,  the  depth  of  the  re- 

f-rcsed  pocket  h„  is  gieetur  than  or  at  least  equal  to  h.  Therefore,  it  is  a  good 
K. 

approximation  to  assume  that  the  pressure  in  the  pocket  is  uniform  so  long  as 
d^/dj.  is  not  too  small. 

f  K 


Now,  we  again  apply  Equation  (2.11)  to  the  inlet  to  the  bearing  film  (see  Figure  1) 
Pn'  V  3h 

pf-Pr-kJV~  (2-13) 


In  solving  the  Reynolds'  equation  (2.4)  with  small  periodic  variations  of  the 
gap  about  the  equilibrium  position,  we  write  in  complex  form, 

H  ■  1  +  e  elT  (2.14) 


and  expand  the  dimensionless  pressure, 

-  .  -  -  It 

p  -  Po  +  e  Pi  e 


12.15) 


taking,  of  course,  only  the  real  part  to  be  of  physical  significance. 

order  in  e  we  can  compute  easily, 

“\ 

P1  “  PQ2  +  2 e  pQ  Pi  ei 
ph  -  +  e  (pQ  +  pi)  e1’ 

h3  "  1  +  3c  e^1 


rp  - 

4. 


l  x  r  i. 


(2.16) 


Substitution  of  the  above  equations  into  (2.4)  yields  the  steady  state  and  perturba¬ 
tion  equations. 


1 

r  Sr 


0 


(2.17) 
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-  —  r  ~  (p  Pl)  -  i  o  (p  +  pt) 
^  3?  I  3r  0  J  0 


(2.18) 


The  boundary  conditions  are; 


at  r  »  1,  pQ  -  1,  Pi  -  0 


(2.19) 


at  r  *»  r 


R  2R’ 


pO  “  PR0 


(2.20) 


It  is  to  be  noted  that  p  .  is  an  unknown  quantity  to  be  determined.  Further- 

RU 

more,  one  additional  boundary  condition. p;  at  rc  is  required  to  solve  the  problem. 

f  *  K 

Now  we  can  obtain  the  steady-state  solution  of  (2.17)  in  terms  of  p  .  Thus, 

.  Ku 


1  + 

in  rr 


1  _  i 

—  i.n  r 


The  radial  velocity  in  the  film  is,  from  the  Stoke's  equation, 


1  3p  ,  a  h2 

U  ■  T-7^  (2^  -  7—  ) 

2p  or  4 


(2.22) 


Here  z  is  the  axial  coordinate  (normal  to  the  film)  measured  from  the  center  oi  the 
film;  the  bearing  surfaces  are  at  z  «  +  h/2  respectively.  Knowing  the  radial 

velocity  distribution  across  the  film,  we  can  calculate  the  mass  flow  rate  per  unit 
circumferential  length  at  the  film  inlet, 

h/2  PD  l  -  |  h/2  2 

G„  *  /  c„  u  dz  »  —  —  I  (z2  — ?-)  dz 

“  -h/2  *  or  '  rn  -h/2 


-h3  3p2 

24i,Rl  3r 


(2.23) 


Define  a  mean  velocity  V  at  the  film  inlet, 

K 


GR  °  PR  VR  h  “  RT  VR  h 


(2.24) 


Equating  (2.23)  and  (2.24)  results  in 


PR  VR 


»  _  i£  ! 

24p  Ir  | 


h2  pa2  3p2  | 

24y  R  -  : 


(2.25) 
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Recognicing  that 


Pr  "  Pro  +  e  PR1  e 


It 


and 


iii 

3r 


rR  ,r 


+  2t  e 


IT  ><*0  Pj) 


3r 


Equation  (2.25)  becomes 
V, 


C2p 


a  1 


R  24mR 


E  p 


RO 


(l  +  teiT 

Pri  3(p0  p  ) 

2 - +  2E  - | - 

3r 

} 

pro 

rJ 

where 


rap  2 

ro 

Dr 

^2 

1  u 

(2.26) 


(2.2?) 


Now  that  we  have  an  expression  for  VR,  let  ua  apply  the  law  of  mass  conservation 
to  the  recessed  pocket. 


1 


^  I  c\  VR  rR  h  "  PF  VF  rF  (h  +  V 

L 

+  k[PTV  <rR2  -  rFZ)(h-  VJ 


(2.28) 


We  nave  assumed  that  the  pressure  and  hence  the  density  in  the  recessed  pocket  be 
uniform  oven  under  dynamic  condition.  This  requires  that  o  <<  1,  as  shown  in 
Appendix  A. 


Alter  some  algebraic  manipulation.  Equation  (2.28)  is  reduced  to 


c?pa 

rR  1 

u 

2upR 

‘F  fro  11  * 

V  E  V 

Pfi 

1 

5r 

V  -  V 

f 

—  *.  c 

PFC 

1  +  hR 

?R 

fi  +  c  eiT 

l  l 


2  +  2E 


3(Pq  V 


or 


Pbo  +  a  +  V  Plil  E 


K  Fbi  I 


(2.29) 


We  have  used  Equation  (2.26)  and  the  ob',.-  -s  relation 

p"f  "  Pro  +  e  PF1  elT 


(2.30) 


Knowing  V  and  V„  from  Eq.'s  (2.26)  and  (2.29)  respectively,  we  are  now  in  a 
K  F 

position  to  utilize  Eq.'s  (2.12)  and  (2.13).  Substituting  (2.26)  into  Eq .  (2.13) 
and  collecting  terms  of  the  same  power  of  e,  we  obtain 


1_  .-3 


PF0  PR0  “  K  Ql  _  E 


PF1  -  PR1  *  (pF0  -  PR0}  7  -  6  —  +  6E 


-1} 


5r  | 


(2.31) 


(2.32) 


Similarly,  from  (2.12)  and  (2.29), 


».  -  ho  •  -  K  1  m-  e-3 

PF0 


-  PF1  "  (P8  -  PF0) 


(  1  +  3: 

6  +  - r 

L  1  +  h 


3ER  ,  .r  3(p0  pl> 
- +  6E  - 


3r  r. 


(2.33) 


pvi  .  V  -  v  r. 


-^-310  - - —  |  PFQ  +  (1  +  ly  Pn"  E  )  (2.34) 

Prn  rD  L  ✓ 


where 


Op  - 

ra 

2  x  243  n4  (  "r'D^R3 

i*?y  -  -2 

—  (1  +  h^)  ^  «■  geometry  parameter 

\r„ 


(2.35) 


(2.36) 


Steady-State  Solution 


Combining  (2.21)  and  (2.27)  results  in 

, . ‘Vjb 

f.n2  -  ‘ 


Thus,  we  rewrite  (2.31)  and  (2.33)  in  the  form 


PR0  (pF0  ~ 


— z -  <PR02  -  1>3-  o 

a  3  d  4  r  RC 

i 


(2.37) 
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PF0  (ps  "  PF0^  3-4  (pR02  "  ° 

Ai  P» 


(2.38) 


whore 


»  modified  feeding  parameter 


12 


U  /  RT  rp  (C  +  h^)  (-  £n  rR) 

“cV - 


16  u  /~RT 
Kp#  (C  +  hR)" 


1/3 


(2. 39) 


Thus,  it  is  revealed  from  the  analysis  that  an  inherently  compensated  hydrostatic 
thrust  gas  bearing  has  two  controlling  parameters,  namely,  a  geometry  parameter,  F, 
and  a  modified  feeding  parameter^  A^,  which  is  defined  in  (2.39).  The  modification 
is  through  the  factor  ^  88  a  result  u8in8  Vohr'e  experimental 

correlation.  From  the  input  data  For  a  given  configuration,  A.  and  T  can  be 
calculated.  Then,  Equations  (2.37)  and  (2.38)  are  to  be  solved  for  pRQ  and  pRQ. 
Since  the  equations  are  non-linear,  the  computations  using  Hewton-Raphson  method 
(Ref.  9)  are  programmed  on  a  computer.  Having  solved  pR(J  and  pRQ,  we  know  the 
steady-atate  pressure  in  the  pocket  (pR^) ,  and  the  pressure  distribution  in  the 
film  is  given  by  Eq .  (2.21). 


From  previous  experience  (see,  for  example.  References  1  and  2),  we  know  that, 
in  order  for  a  bearing  to  be  stable  and  to  have  relatively  large  load  capacity  and 
stiffness,  the  recessed  pocket  should  have  a  relatively  small  volume  but  large  area 
(■rr  2) .  Thus,  a  shallow  pocket  Ol  >  C)  is  desirable.  Under  these  conditions,  it  is 
found  that  T  is  of  the  order  of  1000  or  largar.  From  Equation  (2.37),  it  is  seen  that 
when  F  is  large  (say,  100  or  larger),  pRQ  £  pRQ,  i.e.,  the  loss  at  the  entrance  to  the 
bearing  is  negligibly  small.  Therefore,  as  long  as  F  is  large,  the  solution  will  be 
insensitive  to  F. 


Perturbation  Solution 

The  perturbation  pressure  is  governed  by  Eq.  (2.18)  with  the  boundary  condition 
that  p;  *  !)  at  r  ■  1.  Equations  (2.32)  and  (2.34)  are  to  be  used  to  derive  one 
additional  boundary  condition  required  to  solve  the  problem. 
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Since  pi  may  be  complex,  it  is  convenient  to  assume  that 
PQ  Pi  “  u  +  iv 

Substitute  into  Eq.  (18)  and  separate  the  real  and  imaginary  parts 


(2.40) 


(2.41) 


(2.42) 


with  boundary  conditions 

-  0 


r-1 

-  0 

r-1 

1 

2  - 

PF0  j 

?R 

\ 

PR0  / 

L 

r  dr 


_  t  j  +  F  „ 

-  1  d?  -  r  r 

rR  rR 


+  Fi  Ni  ~  7  PRO  ^pFO  “  PRO^ 


L  PF0  ^ 

1  T  i 

*  Is  - 

+  Li~ 

?"?R 

1  PR0  i 

r  dr  ; 

1  dr 
rR 

+  F  N.  -  F,  N 
r  i  i  r 


(2.43) 


The  last  two  boundary  conditions  are  derived  from  Eq.'s  (2.32)  and  (2.34).  The 

details  of  derivation  and  the  definitions  of  L  ,  L . ,  F  ,  etc.,  are  shown  in 

r  i  r 

Appendix  B. 


In  appendix  C,  Eq.'s  (2.41)  and  (2.42)  subject  to  boundary  conditions  (2. 43.) 
are  solved  numerically  using  matrix  multiplication  method  (Ref.  8). 
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3 .  LOAD  CAPACITY  AND  DYNAMIC  BEARING  REACTIONS 

The  pressure  distribution  in  the  thrust  bearing  under  consideration  can  be 
summarized  as  follows : 

i)  The  pressure  in  the  feeding  hole  region  is  uniform  and  steady, 

ii)  The  pressure  inside  the  recessed  pocket  is  uniform  (approxi¬ 
mately)  but  time-dependent, 

iii)  The  pressure  distribution  in  the  film  is 

p  (r,  t)  -  pQ  (r)  +  e  e  — 3 - 

P- 


“  P  (r)  +  e 


u  (r)  cos  t  -  v  (r)  sin  x 

P0(r) 


The  bearing  force  may  be  obtained  by  Integrating  the  pressure  relative  to  the 
ambient,  throughout  the  film.  Thus, 

R 

W  ■  /  (p  -  p  )  2  rtr  dr 
o 

R 

-  trr  2  (p  -  p  )  +  it (r2  -  r  2)  (p  -  p  )  +  2tt  /  (p  -  p  )  r  dr  (3.2) 


Non-dimensionalizing  the  load  by  ,rR2pa,  we  have 

w  “  7^2  "  V  +  (V  -  V)(PFo  -  1  +  ^  Pfi  eiT) 

1  . 

+2  /  (po  -l+ep  eT)rdr 


Steady-State  Load  Capacity  and  Stiffness 

In  Eq.  (3.3)  the  t  ir^-independent  part  alone  contributes  to  the  steady-state 
load  capacity, 


i&T  *  rv2  (PE  “  +  (?r2  -  rF2)(pF0  “  +  1  I  (Pn  "  1)  r  dr 


wR^p  *F  vrs 

*a 


(3.4) 
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Thls  can  be  easily  calculated  with  the  aid  of  Eq.  (2.21)  and  the  solutions, 

p„n  and  port,  of  Equations  (2.37)  and  (2.38).  From  the  load  capacity,  the  static 
r  U  KO 

stiffness  can  be  obtained  by 


Ck 

o 

ffR^p 

*a 


c 

2  AC 


0.5) 


where  the  superscripts  (+)  and  (-)  refer  to  load  capacities  at  C  +  AC  and  C  -  AC 
respectively.  AC  should  be  sufficiently  small-  a  suitable  value  for  AC  is  0.01C. 


Recall  that  ir.  calculating  the  steady-state  pressures,  we  hove  two  parameters, 
namely,  the  geometry  parameter  T  and  the  modified  feeding  parameter  A^. 

Ck 

In  Figure  3,  the  dimensionless  static  jstiffness  ^3-—  ia  plotted  against 

A.  for  a  bearing  with  ■—  m  .002,  -r^1  -  0.5.  —  «*  2  and  T  ■  i.736  x  106.  It  seems 
1  K  K  C 

that  for  p^  *  2,  3  and  A,  the  respective  stiffness  has  a  maximum  value  when  A^  is 
approximately  0.9.  If  we  change  the  geometry  to  make  F  ■  1  x  103,  numerical 
computation  shows  that  the  dimensionless  stiffness  falls  fairly  closely  with  the  re¬ 
spective  curves  in  Figure  3.  This  confirms  the  conclusion  in  the  previous  section 
that  as  long  as  f  is  large,  the  results  should  be  insensitive  to  T. 


Dynamic  Bearing  Reactions 

The  dynamic  bearing  reaction  due  to  axial  vibration  is,  from  the  time-depend¬ 
ent  part  of  Eq.  (3.3) 


hr  ”  c  Re  \  (*r2  “  rp2)  pFi  e  T  +  2  /  Pi  r  dr 


e  Re  |eiT(Uz  +  i 


(3.6) 


where 


-  -  <;R!  -  V>  (fpih  -  2  !  r-  f  <15 


P» 


Dynamic  Stiffness 


(3.7) 
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V 

z 


Dynamic  Damping 


-  Y2 


"H^i"2/  ^r~Tdl 

;  p„ 


(3.8) 


Hare  u  and  v,  the  solutions  of  (2.41),  (2.42)  and  (2.43),  are  solved 
numerically  by  the  matrix -multiplication  method  as  shown  in  Appendix  C.  Thus, 
the  integrals  in  (3.7)  and  (3.8)  can  be  calculated  numerically.  Using  these 
results,  the  dynamic  stiffness  and  damping  are  plotted  against  frequency  in 
Figures  4  and  5  for  different  C.  It  is  seen  that  when  the  frequency  is  low  (u  <  D 
the  dynamic  stiffness  approaches  asymptotically  to  the  value  of  the  static  stiffness 
and  the  dynamic  damping  approaches  zero,  as  can  be  anticipated.  When  the  frequency 
Increases,  Che  dynamic  damping  first  decreases  and  reaches  a  minimum,  then  it  starts 
to  increase  as  shown  in  Figure  5.  The  frequency  at  which  -  0,  is  called  the 
critical  frequency.  These  will  be  useful  in  the  stability  analysis  in  the  next 
section. 
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4.  STABILITY 


In  the  previous  section,  we  have  calculated  the  dynamic  bearing  reactions 
corresponding  to  small  axial  vibrations  about  the  equilibrium  (statically)  position. 
These  information  are  directly  useful  in  determining  the  bearing  stability. 

In  Reference  6,  a  stability  analysis  for  either  a  single  or  two  degree-of-f reedom 
system  was  performed.  The  results  for  a  single  degree-of-freedom  system  are  directly 
applicable;  they  may  be  stated  as  follows: 


Let  be  the  frequency  of  vibration  at  which 


V 

z 


V 

o 


0 


(4.1) 


This  is  the  state  of  neutral  stability. 

,2 


M 


p  ttR^ 

ra 

Cv 


U 


Then,  the  critical  mass  is  given  by 

(4.2) 


A  slight  variation  from  the  state  of  neutral  stability  would  cause  the  system 
to  be  unstable  if  and  only  if 


av 

_ z 

3v 


v 

o 


<5M  >  C 


(4.3) 


av 

where  6M  is  a  small  mass  increment  above  M  .  From  Figure  5,  - —  >  0.  Therefore, 

o  3v  vQ 

in  order  for  the  bearing  to  be  stable,  6M  must  be  less  than  zero  or,  the  bearing 
mass  must  be  kept  below  the  critical  mass. 


Based  on  the  above  and  a  knowledge  of  Uz  and  V  ,  the  values  of  the  critical  mass 
were  calculated  from  Eq.  (4.2);  they  were  shown  in  Figures  6,  7  and  8.  Since  we  are 
dealing  with  bearings  with  subsonic  flow  throughout  the  passage  and  Vohr's  data 
(Ref.  4)  are  essentially  for  low  Mach  number  flows,  we  calculate  the  Mach  number  at 
r  »  Tp  using  Eq .  (2.12)  to  compute  the  velocity.  The  solid  lines  in  Figures  3,  6, 

7  and  8  are  for  Mach  number  M  <  0.9.  Those  segments  of  the  curves  of  Figs.  3,  6,  7 
and  8  shown  in  dashed  lines  correspond  to  M  >  0.9  and  should  be  regarded  with  caution. 

One  of  the  stability  maps  (experimental  data)  of  Ref.  2  is  reproduced  in  Fig.  9 
where  the  critical  depth  of  the  recessed  pocket  is  plotted  against  the  supply  pressure. 
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Tha  present  analysis  predicts  a  critical  depth  of  0.0011  in  at  -  97.5  psia 

which  is  slightly  below  the  experimental  point  in  Figure  9.  This  indicates  that  the 
present  method  yields  concervative  stability  result.  This  is  believed  to  be  in  part 
due  to  the  assumption  of  a  uniform  pressure  in  the  recess,  pocket. 
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CONCLUSIONS 

1.  Vohr's  entrance  restriction  data  can  be  quite  readily  applied  to  analyze 
the  externally  pressurised  thrust  bearing  with  an  inherently  compensated 
restrictor.  This  analysis  uncovers  two  controlling  parameters,  namely, 
the  modified  feeding  parameter  A^  and  the  geometry  parameter  r ;  the  latter 
represents  the  relative  degree  of  restriction  between  the  exit  of  the  feed¬ 
ing  hole  and  the  exit  of  the  recessed  pocket. 

2.  Steady-state  load  capacity  and  stiffness  were  calculated.  It  was  found 
that  the  static  stiffness  has  a  maximum  value  when  the  modified  feeding 
parameter  is  approximately  0.9. 

3.  Applying  the  stability  theory  of  Ref.  6,  stability  maps  were  constructed 
for  a  particular  bearing  geometry  at  different  supply  pressure.  Experi¬ 
mental  results  of  Ref.  2  were  compared  with  the  present  analysis;  it  was 
found  that  the  result  of  the  present  analysis  is  on  the  conservative  side. 
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6.  RECOMMENDATIONS 

1.  Use  the  well-known  nozzle  formula  to  replace  the  loss  coefficient 
formulation,  and  then  carry  out  the  analysis.  Compare  the  results. 

2.  Perform  limit  cycle  analysis  according  to  the  scheme  suggested  in 
Appendix  D. 

3.  Perform  experiments  to  obtain  more  extensive  data  on  loss  coefficient 
at  higher  supply  pressure.  This  will  take  the  compressibility  effects 
into  account  and  thus,  modify  the  results  for  bearings  with  high 
entrance  Mach  number  (but  still  subsonic) . 
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The  left  hand  Mide  of  Eq .  (A-5)  is  essentially  the  mass  flow  rate  through  the  bearing, 
which  is  constant  everywhere  under  quasi-static  condition  <<  1) .  Therefore,  we 
deduce  that 


(A-6) 
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Since  Is  greater  than  or  at  least  equal  to  C,  the  gap  of  the  pocket  la  at  least 
twice  of  the  bearing  film  thickness.  Hence  we  conclude  that,  if  o  <<■  1,  then  the 
gradient  of  p2  is  much  flatter  in  the  pocket  than  that  in  the  bearing  film,  or, 
the  pressure  in  the  recessed  pocket  can  be  assumed  to  be  essentially  uniform 
(spatially) . 


r 


APPENDIX  B 


The  Derivations  of  Boundary  Conditions 
and  Che  Pocket  Perturbation  Pressure 


From  Eq .  (2.34),  we  can  solve  for  p 


1  +  3  h  d(p  pi) 

- -  +  6  +  E  6  - 2_ -  _ 

1  +  h_  dr  rr 


-  3  i  o  p. 


r  2  -  r  2 
rR  rF 


F  +  i  F. 
r  i 


where 


Ps  “  PF0  PF0 


F  »  3  a  E 


r  ^  -  r 
R  F 


a  +  5^ 


Ox ,  we  can  write 


o  "  (d  )  +  i  (d  ) 

PF1  'PF1  r  PF1  i 


if  we  denote 

a  f  +i,  f. 

r  r  r  1  i 

^Fl^r  F  2  +  F  2 

r  i 


(pFl^i 


-1  F.  +  l.  F 
r  i  i  r 

F  2  +  F  2 
r  i 


•  I  l+3h. 

t  +  - -  +  6 

C  d?|-  1  +  fa 

rR  R 


6  E 


r  2  -  r  2 
v  "  P-  F 

J  E  ?F0 

r„ 
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Subatituting  (B.l)  into  (2.32)  and  solving  for  vs  obtain 


,  PF0  “  1  -  -  ,  d  -  J 

2  "  -  PR1  “  F  +  i  F„  ~  PF0  +  PR0  6  E  -  ^Po  Pl^ 


_i _ \lll_ 

+  1  Fi  L  1  +  E 


A  .  ,  „  .  _  -  rR  "  rF 

—  +  6-3  i  o  E  pFQ  — - 

\  rR 


7  (PF0  "  PR0^ 


Multiply  both  aidaa  by  p 


2  -  ?f-\  <PR0  PR1>  -  6  E  -  PF0  +  PR0  7  +  1  7 

1  Pro /  \  r  1  Mdr  dr/  5 


1  (tT^*  ‘)  F1  +  3  ”  E  ;F0  lR  rF  \J-  7  *so 


FrJ  ~  7  PR0  (pF0  “  PR0J 


Tha  product  pRQ  pRI  can  ba  writtan  as 


Pro  pri  ”  (po  Pl)-  ’  (u  +  lv)- 


Now  1st  ua  daflne 

L_  .  I-.-A- 


r  *  >  Fr‘  +  f'7  '  PF0  +  PR0 1  6  E  p; 


'i  *  “  F  2  +  F.2  6  E  PI 
r  i 


1  +  3 


Ni  -  -  3  o  E  ?F0 


Fr  +  Fi 


z  2  -  i  2 
R  F 


v  *  +  v  <■ 
r  i 
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Chan  Eq.  (B.5)  becomes 


FO 


12  -  ~|  (u  +  i  v) 

Pro/ 


“  L  du 

-  L.— 

+  i  (  L  dv 

1  r  ~ 
-  L  dr 

l'« 

_  *■>£.!] 


?  N  +  F.  N,  +  1(F  N.  -  F.  N  ) 
rr  ii  vri  i  r  J 


7  pro  (pfo  ~  Pro*1 


(B.8) 


Separating  the  real  and  imaginary  parts,  we  obtain  easily  the  following  boundary 
conditions 


2  ~  Pf° 

-  L  4* 

-  L.  ^ 

- 

r  dr 

1  dr 

rR 

1  RO 

r2  .ha} 

-  L  & 

rR 

+  L,  — 

pro 

r  dr 

-  1 
rr 

+  zL 


+  z2 


(B.9) 


(B . 10) 


where 

Z1  "  Fr  Nr  +  Fi  Ni  -  7  Pro  <p'fo  ■  Pro> 

*2  -  Fr  Ht  -  F,  Nr 
Pocket  Perturbation  Pressure 


(B.ll) 


The  perturbation  pressure  at  the  pocket  ppp  can  be  calculated  from  Equations 
(B.3)  and  (B.4).  However,  it  is  important  to  note  that,  when  the  frequency  is  very 
low,  F^  £  0  ,  (in  the  limit  of  zero  frequency,  F^  ■  0;  this  corresponds  to  a  static 
perturbation)  and  from  Eq.  (B.4)  we  have 


(pF1)r 


(B.12) 


For  a  given  geometry  and  a  supply  pressure,  by  varying  the  equilibrium  film 
thickness  one  would  obtain  a  steady-state  pocket  pressure  equal  to  one  half  the  supply 
pressure  at  a  particular  value  of  C.  It  is  obvious  from  Equation  (B.2)  that  when 
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1  - 

PF0  "  2  V 

F  -  0  (B.13) 

r 

Equations  (B.12)  and  (B.13)  indicate  that  undar  the  given  supply  pressure, 

(pF^)^  has  a  singular  behavior.  It  can  be  shown,  however,  by  a  static  perturbation 

analysis,  that  whenever  Ff  -  0,  ■  0.  Therefore,  the  singularity  Is  probably  a 

removable  one,  l.e.,  the  ratio  *>r/Fr  will  approach  a  definite  value.  (One  can  argue 

that  physically  (PF^)r  should  approach  a  definite  value).  But,  the  above  computation 

for  (p  .)  cannot  be  programmed  on  a  computer  because  it  is  unlikely  that  the 
Fir 

machine  will  give  us  the  correct  limiting  value. 


-23 


APPENDIX  C 

The  Matrix  Multiplication  Method  in  Solving  Ordinary 
Dlffarandal  Equations  with  ''Two-Point"  Boundary  Conditions* 


Equations 
may  b«  vri ttan 

fu"  + 


D.E. 


v"  + 

at  x 


B.C. 


at  x 


(2.41)  and  (2.42)  togather  with  their 
in  tha  following  general  form: 


*1  u* 

+  f2  U  + 

f3  v"  +  f„  v'  + 

1 

> 

to 

81  v 

+  g2  V  + 

g3  u"  +  g4  u*  + 

g5  u  - 

"  X!  , 

u  “  Cl  V 

l 1  +  C2  V1  +  c3 , 

V  ■  C4 

«  x2. 

u  ■  F, 

v  -  G 

boundary  conditions  (2.43) 


u 

(C.l) 

88 

(C.2) 

u1  +  C5  y'  + 

(C.3) 

(C.4) 

Hare  we  denote  the  independent  variable  by  x.  The  primes  represent  derivatives 
with  respect  to  x.  The  symbols  fj ,  f2,.  .  .  etc.,  are  known  functions  of  x;  cj ,  c2, 
.  .  .  ate.,  are  known  constants. 


In  central  difference  form,  we  can  write 
k 


u(xk) 


u'(xk) 


u  (xk) 


k+1  k-1 

u _ -  u 

24 

k+1  „  k  ^  k-1 

u  -  2  u  +  u 

72 


(C.5) 


We  have  assumed  that  there  are  N  divisions  between  xj  and  X2,  so  that  k  ■  0,  1,  2, 

,  .  x?  -  Xl 
N,  and  4  •  — L 


Now,  Eq.  (A.l)  take  the  form 
,  k 


k+1  i 
U  1* 


77  + 


24 


+  v 


k+1  | f3k  fuk 
F~  +  24 


+  u 


+  V 


-2  ,  ,  k 

F  +  f2 


F 


fii' 

24 


-2f 


+  £5"!  +  V 


k-1 1  f 
4 


k  k 


24 


■  if 


(C.  6) 


*The  method  of  computation  described  in  this  appendix  is  due  to  Castelli  and 
Pirvics  (Ref.  8). 
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A  similar  equation  may  be  obtained  from  Eq.  ( C. 2),  which  together  with  (C.6)  can 
be  written  in  the  following  matrix  form, 


k  k+1  k  k  0k  k-1  ,k 
Ay  +  B  y  +  C  y  •  d 


(C.7) 


where 


.(  ) 


.(  ) 


,<  ) 


^  +  2A 


k  k 

*a-  +  fc*- 

A^  2A 


-2  ,  k 

77  +  *2 


A7  ~  2A 


k  k 

ii _ 4k_ 

A*  2A 


rt  kn 


86 


k  k 
-f4-  +  lit- 
A~  2A 


1  k 
i_  +  &i- 

F  2i 


-2f  .  ,  k 

—77-+  «S 


2gg  .  k  -2  .  i 

+83  77  +  82 


f  k  f  k' 
li-  -  1*l_ 
A*  2A 


A7  "  2l~ 


(C.8) 


Assume  that  the  y-vector  at  station  T'k  +  1  1  can  be  expressed  by 
k  ,  k 

y  +  m 


,k+1  -  ^  k  '  k 


(C.9) 


where  M  is  an  unknown  matrix  and  m,  an  unknown  vector.  From  (C.9)  we  can  write 
formally 
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k  Jt-l  k-1  .  k-1 

y  ■  M  y  +  a 

and 

k-1  Jc-2  k-2  ^  k-2 

y  ■  m  y  +  m 

Substitute  (C.9)  and  (C.10)  into  (C.7), 

,.k  uk  j.  k  .  _k  k-1  .k  .k  k 

(A  M  +  B  )  y  +  C  y  «  d  -Am 

Thus, 

yk  -  [Ak  ^  +  Bk]_1  [-Ck  yk_1  +  (dk  -  Ak  mk) ] 

Comparing  (C.ll)  with  tha  first  aquation  of  (C.10),  wa  find 
M*1"1  -  [Ak  Mk  +  Bk]  _1  [-Ck] 
mk_1  -  [A^  +  Bk]  -1  (dk  -  Ak  mk) 

Using  (C.4) 

F 

N  f 

y  ”  _ 

G 

and  from  the  first  equation  of  (C.10)  wa  obtain 

m"-1  -  0 


Now  we  can  use  (C.12)  and  (C.13)  as  recurrance  formulas  to  obtain 


M*1-2  - 

[a11-1  tf-1  A  A11-1] 

-| 

- 

N-2 

,N-1  WN-1  ,  „N-l[ 

-i  j 

jN-1 

,N-1  N-l 

m  *  1 

-A  M  +  B  J 

L  d 

-  A  m 

m"-3  -  1 
N-3 

m  • 

•A1*-2 

L»-2  h21-2  +  b"-2] 

r 

i 

UH 

[dN-2 

. 

.N-2  N-2| 

-  A  m  J 

- 

(C.10) 


(C.ll) 

(C.12) 

(C.13) 

(C.14) 

(C . 15) 

(C.16) 

(C.17) 

(C.18) 


and  so  on. 
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Having  •- ^reputed  cha  M's  and  n'i,  wa  can  calculato  the  solution  by  marching  from 
x  -  xi  rut  first  1st  us  rewrite  boundary  condition  (C.3). 

u°  ■  Cj  (u1  -  u°)  +  C2  (v1  -  v°)  +  C3 

v°  -  Cl  (u1  -  u°)  +  C5  (v1  -  v°)  +  C6 

using  forward  dlffarenca  formula.  In  matrix  form  tha  abova  two  equations  can  be 
written  as 

S  y°  -  T  y1  +  2  (C.19) 

where 


1 


Using  the  relationship 

y1  ■  M°  y°  +  m°  (C.21) 

Equatio  ,  (C.19)  becomes 

[S]  y°  -  [T][M°]  y°  +  [T]  m°  +  z  (C.22) 

Thus,  we  have  finally 


y°  -  [S  -  x  m°]-1  (T  m°  +  z) 


(C. 23) 
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Now  wc  can  starch  from  this  point  to  obtain  tha  solution,  using  Eq.  (C.9), 


y1  -  M°  y°  +  m° 
y2  -  M1  y1  +  m1 


(C.24) 


and  so  on. 
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APPENDIX  D 

Scheme  for  Limit  Cycle  Analysis 

Let  us  symbolically  write  the  Reynold's  equation  in  the  form: 

R  (f)  -  !s  s’  i  If)  -  (5  s>  '  0  <D1> 


Suppose  that  the  bearing  is  unstable  so  that  the  amplitude  of  oscillation  becomes 
larger  and  larger.  Wow,  let  us  assume  that  a  limit  cycle  motion  will  be  reached 
and  during  the  limit  cycle  motion  (finite  amplitude),  the  gap  variation  and  pres¬ 
sure  fluctuation  remain  to  be  sinusoidal.  Thus, 


h  =  1  +  e  cos  vt  (D.2) 

P  *  PQ  +  P^  (r)  cos  vt  +  qL(r)  sin  vt  (D.3) 


Here,  the  dimensionless  amplitude  c  is  not  necessarily  small. 
Apply  cosine  and  sine  transforms  to  (D.l). 


.*/v 

V  VJ 

>  =  / 

■'-n/v 

R  cos  vt  d  t  =0 

(D  4) 

> 

Tt/  V 

V  v 

J 

1  ■  L 

R^p^sin  vt  dt  =  0 

(D.5) 

Note  that  a  typical  quantity  in  (D.l),  h  ,  is  equal  to 

.3  .  3 

h  *  (1  +  c  cos  vt) 

.22  .33 

=  ]  +  3e  cos  vt  +  3e  cos  vt  +  e  cos  vt 


(D.6) 
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Utilizing 

2 

cos 

3 

cos 


the  identities 

vt  =  1|  (1  +  cos  2vt) 

vt  =  i,  (3  cos  vt  +  cos  3vt) 

cos  mvt  cos  nvt  dt  *=  n  6 

mn 

cos  mvt  sin  nvt  dt  =  0 


sin  mvt  sin  nvt  dt  =  /(  6 

mn 


etc.,  one  realizes  that  (D.4)  and  (D.5)  are  two  non-linear  ordinary  differential 
equations  in  and  q^,  which  are  truncated  to  include  only  the  simple  harmonics. 

Similarly  we  can  carry  out  the  cosine  and  sine  transforms  to  the  boundary  con¬ 
ditions  and  mass  conservation  in  the  recessed  pocket  and  so  on.  This  would  yield 
boundary  conditions  for  the  non-linear  differential  equations  (D.4)  and  (D.5). 
Iteration  method  may  be  used  in  solving  the  non-linear  system. 


Having  solved  p,  end  we  can  compute  the  bearing  forces. 


W  =  2n  J 


W  =  2n  /  (p  -  p  )  r  dr  =  W  +  U  cos  vt  +  V  sin  vt 
a  o 


Then,  the  power  output  of  the  gas  film,  E,  is 


(D.7) 


E  = 


W 


oh 

ot 


oh 

6t 


is 


obv  ious ly 


the 


squeeze  velocity. 


(D.8) 


where 
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The  average 

power  output  over 

a  squeeze  cycle  can  be  readily  obtained  by 

n/v 

n/v 

V 

E  *  ■— 

f  E  dt  *  — 

r 

w 

T1  dt 

in 

1  2n 

j 

dt 

-n/v 

-n/v 

n/v 

»  — 

2n 

(-  cCv)  / 

W  sin 

vt 

dt 

-n/v 

vcC 

2 


V(i,v) 


(D.9) 


Let  the  bearing  mass  be  M  and  write  the  equation  of  motion 


i 

M  — r  +  U  cop  '/t  +  V  sin  vt  ■  F  sin  vt  (D.10) 

dt 

where  F  sin  vt  is  the  externcl  force  acting  on  the  beari-iq  other  than  the  gas  film 
forces . 


d  ^  h  2 

Sinc3  — r  =  -  eCv  cos  vt,  by  collecting  the  cosine  terms  in  (D.10)  we  obtain 
dt 


n 

w 


(L>.  11  j 


Note  that  U  is  a  function  of  v  and  e.  Therefore,  (D.ll)  is  a  highly  non-linear 

equation.  From  equation  (D.ll)  the  natural  frequency,  v  ,  can  be  determined 

c 

using  numerical  methods,  In  general,  v  is  a  function  of  c  and  M.  Thus, 

o 

V  =  V  (e,M)  (D.P) 

o  o 

Having  determined  the  natural  frequency,  a  stability  criterion  i;.,  hr.  established 
by  energy  consideration.  The  bearing  is  stable  if  E  is  negative,  and  unstable  if 
E  is  positive.  From  Equation  (D.9)  we  conclude  that  the  bearing  is  stable  if  V 
is  positive,  and  unstable  if  V  is  negative. 
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Now,  at  v  =  v  , 
o 


V  (€fv) 


-  V 


v 

I  o 


(€,vo(e,M)| 


V  (e,M) 


(D-13) 


A  qualitative  plot  of  V  versus  €  for  a  given  M  is  believed  to  be  that  shown  in 
the  diagram.  The  bearing  is  ut stable  when  e  is  small.  Hence  the  amplitude  will 


grow  until  it  reaches  a  value,  e  ,  where  V  *  0,  and  the  bearing  becomes  morgin- 

o 

ally  stable.  The  bearing  is  then  said  to  reach  a  limit  cycle  motion  with  limit 
cycle  amplitude  e . 
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NOMENCLATURE 

C  equilibrium  film  thickness 

d^.  diameter  of  ^feeding  hole 

dp  diametaf  of  recessed  pocket 

E  defined  in  (2.27)  ^  , 

defined  in>(B.^2) 

Fz  dynamic  bearing  reaction  :  _ 

G  mass  flow  rate  per  unit  length  -  >  ■ 

h  film  thickness  ■  . 

h  -depth  of  recessed  pcfcket  » 

R  ,  ■ 

&  ■  h/C  ’ 

i  ■  •  * 1 

V 

k  static  stiffness 

.  o. 

K1  ,  .  loss  coefficient  -  K  Re 

Lr>Li,  defined  in  (B.7) 

ir  £  -  defined  in  (B.4)  * 

L  ■  typical  length 

M  mass;  matrix  defined  in  (C.9) 

Nr,Ni  defined  in  (B.7) 

p^  pressure 

pfl  supply  pressure 

Pa  ambiant  pressure 

p  dimensionless  pressure 

p  ,pj  defined  in  (2.15) 

O  i 

Ql  defined  in  (2.35)  . 

r  radial  coordinate 

r  dimensionless  radial  coordinate,  r/R 

Tp  radius  of  feeding  hole 

r_  radius  of  recessed  pocket 

R  bearing  radius 

Re  Reynolds'  number  ■ 

Re  Film  entrance  Reynolds'  number  "  — —  -  2 

"V  V 

R  gas  constant 

■ 

T  temperature 

t  time 


dynamic  stiffness 
dynamic  damping 
defined  in  (2.40) 
velocity 
bearing  load 
static  load 

Ul  ,  22]  defined  in  (B.ll) 

geometry  parameter  defined  in  (2.36) 
amplitude  of  axial  vibration 
dimensionless  £ 

modified  feeding  parameter  defined  in  (2.39) 

viscosity 

density 

squeeze  number  defined  in  (2.5) 

dimensionless  time 

frequency 
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